Effects of intrinsic spin-relaxation in molecular magnets on current-induced 

magnetic switching 
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Current-induced magnetic switciiing of a single magnetic molecule attached to two ferromagnetic 
contacts is considered theoretically, with the main emphasis put on the role of intrinsic spin re- 
laxation processes. It is shown that spin-polarized current can switch magnetic moment of the 
molecule, despite of the intrinsic spin relaxation in the molecule. The latter processes increase the 
threshold voltage (current) above which the switching takes place. 

PACS numbers; 75.47.Pq, 75.60.Jk, 71.70.Gm, 75.50.Xx 
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Single-molecule mag nets (SMMs) [H, Q attract much 
attention due to their exceptional properties and possi- 
ble applications in quantum information processing [sj 
and information storage technology [l|, Apart from 
this, SMMs are also promising as key elements of novel 
spintronics devices [5|]. Therefore, an important question 
is how to manipulate the SMM in order to write a bit 
of information on it. One possibility relies on the ap- 
plication of an external magnetic field. In the following 
paper, however, we consider another possibility, i.e. the 
current-induced magnetic switching (CIMS) 0,0]. The 
phenomenon of CIMS 'Sil is well known in the case of 
artificial layered nanostructures. Since the present-day 
technology allows to attach a SMM to electronic con- 
tacts [9] , CIMS of a SMM is an alternative way of writing 
information in SMM-based memory elements. 

There are several challenging aspects of the current- 
induced manipulation of SMM's spin. First, the up-to- 
date experimental techniques offer only limited control 
of the relative orientation of the molecule's easy axis 
and leads' magnetizations [loj . Second, intrinsic spin- 
relaxation time of the molecule [ill has a significant in- 
fiuence on the switching parameters and is hardly con- 
trollable externally. Finally, the efficiency of spin injec- 
tion from ferromagnetic leads to molecules is a subject of 
intense technological efforts. The main objective of this 
paper is a detailed analysis of the second point, i.e. of 
the influence of intrinsic spin-relaxation on the CIMS of 
a SMM. 

It is only very recently, when the switching of SMM's 

Sin due to spin-polarized current has been proposed J6, 
, However, the intrinsic spin relaxation in the molecule 
has not been taken into account. When the energy e 
of the lowest unoccupied orbital (LUMO) level of the 
molecule is sufficiently low, electronic transport takes 
place owing to tunneling between the electrodes and the 
LUMO level. The CIMS can then occur when the LUMO 
level is exchange coupled to the SMM's spin. The corre- 
sponding Hamiltonian of the molecule can be written in 



the form 

HsMM = -{D + J2Di 4 + D2 cjc^cjcj^2 

(T 

+ ^£ciCa + ~ \ ■ '^4co-', (1) 

where cr is the Pauli spin operator for electrons in the 
LUMO level, c\{c„) is the relevant creation (anihilation) 
operator, and U is the Coulomb energy of two electrons 
of opposite spins in the LUMO level. The first term of 
Ti-sMM describes the anisotropy of a SMM, whereas the fi- 
nal one accounts for the exchange interaction between the 
SMM's core and the LUMO level, with J being the rele- 
vant exchange parameter. The infiuence of molecule's ox- 
idation state on the anisotropy [l2| is taken into account 
by the terms linear in Di and D2. In turn, tunneling 
processes between the molecule and leads are described 

] , where Tq 



kcr 



by T^T, Ht = E,Eka [T.atc, 
is the tunneling matrix element between the SMM and 
the (7-th lead {q = L{R) for the left (right) electrode), 
and a^^ {o^ka) annihilation (creation) operator of 

an electron with the wave vector fc and spin a in the 
g-th electrode. The system is shown schematically in 
Fig. Ha). ^ 

Tunneling between the leads and molecule gives rise 
to a finite spin-dependent width of the LUMO level, 
= Eq r« , where T% = 2TT\Tq\^D% and is the spin- 
dependent density of states (DOS) at the Fermi level in 
the lead q. The parameters T% will be used in the fol- 
lowing to describe coupling strength between the LUMO 
level and leads. It is convenient to write as F^ = 
F,(l ± P,), where F, = {T\ + Tl)/2, and P, is the po- 
larization of the q-i\i lead, P, = {D\ - Dl) / {D\ + D'i_) . 
Here cr = +(— ) corresponds to spin-majority (spin- 
minority) electrons. In the following, we assume that 
the couplings are symmetric, Fl = Fr = F/2. 

When the energy e of the LUMO level is large enough, 
electron tunneling to the molecule is energetically forbid- 
den at bias voltages of interest. However, current still can 
fiow due to higher order processes, e.g. cotunneling ones. 
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FIG. 1: (color online) (a) Schematic representation of the sys- 
tem and switching mechanism due to spin-polarized current, 
(b) Energy levels of the Mni2 molecule for the following pa- 
rameters: D « 0.05 mV, Di -0.006 meV, D2 ^ 0.0017 
meV [il, J = 0.25 meV, £ = 5 meV, and U = 0. Different 
parabolas correspond to indicted values of the SMM's total 
spin St and occupation numbers of the LUMO level. 



and CIMS of the molecule's spin is still possible f6| when 
the electrons virtually entering the molecule couple to the 
molecule's spin via the exchange interaction. The Hamil- 
tonian of the molecule can be then reduced to TLsmm = 
—DSl, while tunneling processes can be described effec- 
tively hyUT ^ \ T.qq' T.aa'kk-iJf^'yy' ■ S + '^)«fcl4v" 
where 5 takes into account those tunneling processes be- 
tween the leads, which are not included in the exchange 
term. These, however, are irrelevant from the point of 
view of switching process and can be neglected [5 = 0). 

Switching of the SMM's spin takes place consecutively 
via the magnetic states of the molecule. These states are 
described by the eigenvalue m of the z component of the 
molecule's total spin, S^ = 5'z-|-i(c|c-|-— cjc^) (where the 
second term represents the contribution from electrons 
in the LUMO level), and the corresponding occupation 
number n of the LUMO level, i.e. |n, m) |7|. 

The energy levels corresponding to the Hamiltonian 
Ti-SMM form the four parabolas shown in Fig. [Tl^b). The 
lowest parabola corresponds to the empty LUMO level 
and is the only relevant energy spectrum in the limit 
of high LUMO level. In the latter case, the spin rever- 
sal proceeds via all the consecutive intermediate states 
10,-5"+ 1).. .10,5"- 1) (dots in Fig.IHb)). A more com- 
plicated situation occurs when electrons tunnel directly 
through the LUMO level, and transitions between neigh- 
boring molecular states are governed by the following se- 
lection rules |AS'j^| — 1/2 and the oxidation state of 
the SMM changes by one. Consequently, direct transi- 
tions between molecular states of the same spin multiplet 
are forbidden. 

A convenient way of analyzing the magnetic switching 



of a SMM is to consider the mean value of the z compo- 
nent of the total molecule's spin. 



(2) 



The problem is then reduced to determining the probabil- 
ities P\n,m) of finding the molecule in all possible molec- 
ular states |n, m). These probabilities can be determined 
from the relevant master equations and the correspond- 
ing transition rates between the molecular states. The 
key point is that these transition rates must include also 
the effects due to intrinsic spin relaxation. 

Generally, in the systems under consideration one can 
distinguish two classes of SMM's spin relaxation pro- 
cesses. The first class is associated with the coupling 
of the molecule to ferromagnetic leads d, 0, [11]) and 
the other one includes all intrinsic spin-relaxation pro- 
cesses 0, [ll|. The role of the latter processes in the 
CIMS of the SMM's spin is the main objective of this 
paper. It is important to note that even at low tempera- 
tures the molecule's spin is subject to decoherence due to 
interaction with its environment. A SMM in an excited 
molecular spin level can undergo transitions to neigh- 
boring levels of lower energy, which is accompanied by 
emission of a phonon. As a consequence, excited molec- 
ular spin states have a finite life-time, and it has been 
shown that this time for Peg is of order of 10~^ s 
Furthermore, coherence of the SMM's spin can also be 
lost due to various forms of magnetic interactions with 
the environment, e.g. due to the hyperfine interaction 
with nuclear moments of protons in the vicinity of the 
molecule P,'Tl'|. 

To include the intrinsic spin relaxation processes into 
considerations, we introduce the relaxation rate in 
addition to the rates 'j'!"-™)!" ) describing current- 
induced transitions between the molecular spin states 
\n,m) and \n',m'). The latter ones can be calculated 
from the Fermi golden rule 0| • In turn, intrinsic re- 
laxation of the molecule's spin occurs as transitions be- 
tween neighboring molecular states of the same spin mul- 
tiplet. Fig. [T]Jb), i.e. the occupation of the LUMO level 
is not changed by these processes. Furthermore, we as- 
sume that the spin relaxation is fully characterized by a 
phenomenological relaxation time r^, i.e. the relaxation 
rate takes the form 
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Here, €\n,m) denotes the energy of the molecular state 
71, m), /cb is the Boltzmann constant, T is the tempera- 
ture of the system, and A = e\n,m) ~ £|n,m')- The Boltz- 
mann factor in Eq. fSj assures that the intrinsic spin re- 
laxation drives the SMM's spin to the state of the lowest 
energy. 
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FIG. 2: (color online) The eflFect of intrinsic relaxation pro- 
cesses on magnetic switching of the molecule Mni2 in the limit 
of high LUMO level, calculated for indicated values of the re- 
laxation time, Tu, and for parallel magnetic configuration. 
The polarization parameters of the electrodes are: Pl = 1 
and Pr = 0.5. The other parameters are: J = 100 meV, 
D « 0.05 meV [H, T = 0.01 K, and c = 10 kV/s. 



Taking into account the relaxation processes discussed 
above, the master equations for the probabilities P\n,m) 
take the form, 



dB 
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In the following we assume that initially the molecule is 
saturated in the state |0, —10), and then voltage growing 
linearly in time is applied, V — ct, with c denoting the 
speed at which the voltage is augmented. It means that 
for the molecule of the spin S = 10, like the molecule 
Mni2 or Feg, one has to solve the set of 21 coupled dif- 
ferential equations for the situation of large LUMO level 
and 84 equations in the general case. 

In Fig. [5] we show evolution of the z component of the 
molecule's spin in the case of parallel magnetic configu- 
ration and high LUMO level (current flows then due to 
higher order processes). The results clearly show that 
the molecule's spin becomes switched when the voltage 
exceeds some critical value, which is determined by the 
magnetic anisotropy (energy gap between the states cor- 
responding to m = — 10 and m — —9) and the intrin- 
sic relaxation time. Since the intrinsic spin-flip relax- 
ation processes tend to restore the initial state, the low- 
est threshold voltage occurs in the absence of intrinsic 
spin relaxation. The switching, however, takes also place 
in the presence of intrinsic spin relaxation processes, al- 
though the threshold voltage becomes increased. Apart 
from this, the switching time also increases with decreas- 
ing Tjj. Similar behavior also occurs in the case when 
magnetic moments of the leads are antiparallel. 

The parameters assumed in Fig. [2] correspond to half- 
metallic ferromagnetic left electrode {P^ — 1), and typ- 
ical 3d ferromagnetic metallic right electrode. For sim- 
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FIG. 3: (color online) (a) The mean value of the total spin 
(St) for different values of the inverse relaxation time, , in 
the case of parallel configuration of the electrodes' magnetic 
moments, calculated for Pl = Pr = 0.5. Solid lines in the 
part (b) represent cross sections of the plot (a) for several 
values of tr, whereas the dashed lines correspond to the cur- 
rent fiowing through the system. The other parameters are 
as in Fig. 1, and T = 0.01 K, c — 1 V/s, and the coupling 
parameter F — 0.001 meV. 



plicity the positive bias corresponds to electrons flowing 
from left to right (e > 0), i.e. from half- metallic fer- 
romagnetic electrode to the 3d one. Spin-up electrons 
leaving the half-metallic electrode can change its spin 
orientation when interacting via exchange coupling with 
the molecule's spin, and this way can increase the spin 
number m of the molecule's spin. Intrinsic relaxation 
processes tend to restore the initial state. When the cur- 
rent exceeds some critical value, the competition of in- 
trinsic spin relaxation (lowering the quantum number m) 
and current-induced processes (increasing the number m) 
leads to spin reversal of the molecule. This takes place 
in both, parallel and antiparallel (with magnetic moment 
of the right electrode being reversed) magnetic configura- 
tions. For reversed bias polarization only switching from 
the state |0, 10) to the state |0, —10) is possible. 

In Figs [3] and [4] we show the average value of the to- 
tal spin (S^) and current flowing in a biased system in 
the case when switching occurs due to sequential tun- 
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FIG. 4: (color online) The same as in Fig. 3 but for antipar- 
allel alignment of the electrodes' magnetic moments. 

neling of electrons through the molecule's LUMO level. 
These two figures correspond to parallel (Fig. [3]) and an- 
tiparallel (Fig. [5]) magnetic configurations. Clearly, there 
is no switching in the parallel configuration. Instead of 
this, current excites the molecule to higher states and 
the average spin becomes zero (see Fig. [3]) . The situation 
is different in the antiparallel configuration, where there 
is a clear switching from the state |0,— 10) to the state 
|0, 10). To understand this behavior one should note that 
in Figs [3] and H] the spin polarization of both electrodes 
is the same. Consequently, the current-induced processes 
increasing the number m and those decreasing m occur 
with the same rate in the parallel configuration. Accord- 
ingly, none of the molecule's spin states is stabilized by 
the current. In contrast, in the antiparallel configura- 
tion processes increasing the number m start to domi- 
nate over those decreasing m above a certain threshold 
voltage, and the switching to the state |0, 10) takes place. 
Current-induced switching of the molecules's spin may be 
possible also in the parallel configuration, provided spin 
polarizations of the electrodes are different. 

In conclusion, we have shown that spin-polarized cur- 
rent flowing though the molecule can switch its magnetic 
moment despite of intrinsic spin relaxation processes in 
the molecule. The latter processes increase the thresh- 



old voltage (current) and switching time. If for a cer- 
tain bias polarization current stabilizes the state |0, —10) 
(or |0, 10)), then the opposite current stabilizes the state 
|0,10) (or |0,-10)). 
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